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Abstract I discuss the use of Type la supernovae (SNe la) for cosmological dis- 
tance determinations. Low-redshift SNe la [z ^ 0.1) demonstrate that 
the Hubble expansion is linear with Ho = 72 ± 8 km s _1 Mpc -1 , and 
that the properties of dust in other galaxies are generally similar to those 
of dust in the Milky Way. The light curves of high-redshift (z — 0.3- 
1) SNe la are stretched in a manner consistent with the expansion of 
space; similarly, their spectra exhibit slower temporal evolution (by a 
factor of 1 + z) than those of nearby SNe la. The measured luminosity 
distances of SNe la as a function of redshift have shown that the ex- 
pansion of the Universe is currently accelerating, probably due to the 
presence of repulsive dark energy such as Einstein's cosmological con- 
stant (A). From about 200 SNe la, we find that H t = 0.96 ± 0.04, 
and Qa — IAVIm = 0.35 ± 0.14. Combining our data with the results 
of large-scale structure surveys, we find a best fit for Qm and Aa of 
0.28 and 0.72, respectively — essentially identical to the recent WMAP 
results (and having comparable precision). The sum of the densities, 
~ 1.0, agrees with extensive measurements of the cosmic microwave 
background radiation, including WMAP, and coincides with the value 
predicted by most inflationary models for the early Universe: the Uni- 
verse is flat on large scales. A number of possible systematic effects 
(dust, supernova evolution) thus far do not seem to eliminate the need 
for Q,a > 0. However, during the past few years some very peculiar 
low-redshift SNe la have been discovered, and we must be mindful of 
possible systematic effects if such objects are more abundant at high 
redshifts. Recently, analyses of SNe la at z = 1.0-1.7 provide further 
support for current acceleration, and give tentative evidence for an early 
epoch of deceleration. The dynamical age of the Universe is estimated 
to be 13.1 ± 1.5 Gyr, consistent with the ages of globular star clusters 
and with the WMAP result of 13.7 ± 0.2 Gyr. According to the most 
recent data sets, the SN la rate at z > 1 is several times greater than 
that at low redshifts, presumably because of higher star formation rates 
long ago. Moreover, the typical delay time from progenitor star forma- 
tion to SN la explosion appears to be substantial, ~ 3 Gyr. Current 
projects include the search for additional SNe la at z > 1 to confirm 
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the early deceleration, and the measurement of a few hundred SNe la at 
z — 0.2-0.8 to more accurately determine the equation-of-state param- 
eter of the dark energy, w = P/(pc 2 ), whose value is now constrained 
by SNe la to be in the range —1.48 £ w — 0.72 at 95% confidence. 

Keywords: cosmological parameters, dark energy, distance scale, supernovae 

Introduction 

Supernovae (SNe) come in two main observational varieties (see Filip- 
penko 1997b for a review). Those whose optical spectra exhibit hydrogen 
are classified as Type II, while hydrogen-deficient SNe are designated 
Type I. SNe I are further subdivided according to the detailed appear- 
ance of the early-time spectrum: SNe la are characterized by strong 
absorption near 6150 A now attributed to Si II, SNe lb lack this feature 
but instead show prominent He I lines, and SNe Ic have neither the Si II 
nor the He I lines (at least not strong ones). SNe la are believed to re- 
sult from the thermonuclear disruption of carbon-oxygen white dwarfs, 
while SNe II come from core collapse in massive supergiant stars. The 
latter mechanism probably produces most SNe Ib/Ic as well, but the 
progenitor stars previously lost their outer layers of hydrogen or even 
helium, through either winds or mass transfer onto a companion star. 

It has long been recognized that SNe la may be very useful distance 
indicators for a number of reasons; see Branch & Tammann (1992), 
Branch (1998), and references therein. (1) They are exceedingly lumi- 
nous, with peak Mb averaging —19.0 mag if Hq = 72 km s _1 Mpc -1 . (2) 
"Normal" SNe la have small dispersion among their peak absolute mag- 
nitudes (<t 0.3 mag). (3) Our understanding of the progenitors and 
explosion mechanism of SNe la is on a reasonably firm physical basis. 
The results of recent models give good fits to the observed spectra and 
light curves, providing confidence that we are not far off the mark. (4) 
Little cosmic evolution is expected in the peak luminosities of SNe la, 
and it can be modeled. (5) One can perform local tests of various possi- 
ble complications and evolutionary effects by comparing nearby SNe la 
in different environments (elliptical galaxies, bulges and disks of spirals, 
galaxies having different metallicities, etc.). 

Research on SNe la in the 1990s has demonstrated their enormous 
potential as distance indicators. Although there are subtle effects that 
must indeed be taken into account, it appears that SNe la provide among 
the most accurate values of Hq, qo (the deceleration parameter), Qm (the 
matter density), and Qa [the cosmological constant, Ac 2 /(3Hq)]. 

For more than a decade there have been two major teams involved 
in the systematic investigation of high-redshift SNe la for cosmolog- 
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ical purposes. The "Supernova Cosmology Project" (SCP) is led by 
Saul Perlmutter of the Lawrence Berkeley Laboratory, while the "High- 
Z supernova search Team" (HZT) is led by Brian Schmidt of the Mt. 
Stromlo and Siding Springs Observatories. I have been privileged to 
work with both of these teams (see Filippenko 2001 for a personal ac- 
count), but my primary allegiance is now with the HZT. A few years 
ago, the HZT split into two overlapping subsets: the "Higher-Z Super- 
nova Search Team" led by Adam Riess of the Space Telescope Science 
Institute, and the ESSENCE team ("Equation of State: SupErNovae 
trace Cosmic Expansion") led by Christopher Stubbs of Harvard Uni- 
versity. Other groups have recently formed to conduct similar stud- 
ies, such as the supernova team of the Canada-France-Hawaii Telescope 
Legacy Survey. An outgrowth of the SCP, the very large SNAP (Su- 
per Nova/ Acceleration Probe) collaboration is planning a future space 
satellite that will be largely devoted to using SNe la for cosmology. 

1. Homogeneity and Heterogeneity 

Until the mid-1990s, the traditional way in which SNe la were used for 
cosmological distance determinations was to assume that they are per- 
fect "standard candles" and to compare their observed peak brightness 
with that of SNe la in galaxies whose distances had been independently 
determined (e.g., with Cepheid variables). The rationale was that SNe la 
exhibit relatively little scatter in their peak blue luminosity (ctb ~ 0.4- 
0.5 mag; Branch & Miller 1993), and even less if "peculiar" or highly 
reddened objects were eliminated from consideration by using a color 
cut. Moreover, the optical spectra of SNe la are usually rather homo- 
geneous, if care is taken to compare objects at similar times relative to 
maximum brightness (Riess et al. 1997, and references therein). Over 
80% of all SNe la discovered through the early 1990s were "normal" 
(Branch, Fisher, & Nugent 1993). 

From a Hubble diagram constructed with unreddened, moderately 
distant SNe la (z ^ 0.1) for which peculiar motions are small and rela- 
tive distances (given by ratios of redshifts) are accurate, Vaughan et al. 
(1995) find that 

(M B (max)} = (-19.74 ±0.06) +5 log (H /50) mag. (1) 

In a series of papers, Sandage et al. (1996) and Saha et al. (1997) 
combine similar relations with Hubble Space Telescope (HST) Cepheid 
distances to the host galaxies of seven SNe la to derive Hq = 57 ± 4 km 
s" 1 Mpc" 1 . 
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Over the past two decades it has become clear, however, that SNe la 
do not constitute a perfectly homogeneous subclass (e.g., Filippenko 
1997a, b). In retrospect this should have been obvious: the Hubble dia- 
gram for SNe la exhibits scatter larger than the photometric errors, the 
dispersion actually rises when reddening corrections are applied (under 
the assumption that all SNe la have uniform, very blue intrinsic col- 
ors at maximum; van den Bergh & Pazder 1992; Sandage & Tammann 
1993), and there are some significant outliers whose anomalous magni- 
tudes cannot be explained by extinction alone. 

Spectroscopic and photometric peculiarities have been noted with in- 
creasing frequency in well-observed SNe la. A striking case is SN 1991T; 
its pre-maximum spectrum did not exhibit Si II or Ca II absorption lines, 
yet two months past maximum brightness the spectrum was nearly in- 
distinguishable from that of a classical SN la (Filippenko et al. 1992b; 
Ruiz-Lapuente et al. 1992; Phillips et al. 1992). The light curves of 
SN 1991T were slightly broader than the SN la template curves, and 
the object was probably somewhat more luminous than average at max- 
imum. Another well-observed, peculiar SNe la is SN 1991bg (Filippenko 
et al. 1992a; Leibundgut et al. 1993; Turatto et al. 1996). At maximum 
brightness it was subluminous by 1.6 mag in V and 2.5 mag in B, its 
colors were intrinsically red, and its spectrum was peculiar (with a deep 
absorption trough due to Ti II). Moreover, the decline from maximum 
was very steep, the /-band light curve did not exhibit a secondary max- 
imum like normal SNe la, and the velocity of the ejecta was unusually 
low. The photometric heterogeneity among SNe la is well demonstrated 
by Suntzeff (1996) with objects having excellent BVRI light curves. 

2. Cosmological Uses: Low Redshifts 

Although SNe la can no longer be considered perfect "standard can- 
dles," they are still exceptionally useful for cosmological distance deter- 
minations. Excluding those of low luminosity (which are hard to find, 
especially at large distances), most of the nearby SNe la that had been 
discovered through the early 1990s were nearly standard (Branch et al. 
1993; but see Li et al. 2001b for more recent evidence of a higher intrinsic 
peculiarity rate). Also, after many tenuous suggestions (e.g., Pskovskii 
1977, 1984; Branch 1981), Phillips (1993) found convincing evidence for 
a correlation between light-curve shape and the luminosity at maximum 
brightness by quantifying the photometric differences among a set of 
nine well-observed SNe la, using a parameter [Amis(-B)] that measures 
the total drop (in B magnitudes) from i?-band maximum to t = 15 days 
later. In all cases the host galaxies of his SNe la have accurate relative 
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distances from surface brightness fluctuations or from the Tully-Fisher 
relation. The intrinsically bright SNe la clearly decline more slowly than 
dim ones, but the correlation is stronger in B than in V or /. 

Using SNe la discovered during the Calan/Tololo survey (z 0.1), 
Hamuy et al. (1995, 1996b) refine the Phillips (1993) correlation between 
peak luminosity and Ami^(B). Apparently the slope is steep only at 
low luminosities; thus, objects such as SN 1991bg skew the slope of the 
best-fitting single straight line. Hamuy et al. reduce the scatter in the 
Hubble diagram of normal, unreddened SNe la to only 0.17 mag in B 
and 0.14 mag in V; see also Tripp (1997). Yet another parameterization 
is the "stretch" method of Perlmutter et al. (1997) and Goldhaber et al. 
(2001): the B-band light curves of SNe la appear nearly identical when 
expanded or contracted temporally by a factor (1 + s), where the value 
of s varies among objects. In a similar but distinct effort, Riess, Press, 
& Kirshner (1995) show that the luminosity of SNe la correlates with 
the detailed shape of the overall light curve. 

By using light-curve shapes measured through several different filters, 
Riess, Press, & Kirshner (1996a) extend their analysis and objectively 
eliminate the effects of interstellar extinction (as is also now done with 
the Ami5 method; Phillips et al. 1999). A SN la that has an unusually 
red B — V color at maximum brightness is assumed to be intrinsically 
subluminous if its light curves rise and decline quickly, or of normal lumi- 
nosity but significantly reddened if its light curves rise and decline more 
slowly. With a set of 20 SNe la from the Calan/Tololo sample (Hamuy et 
al. 1996c) and the Harvard/Smithsonian Center for Astrophysics (CfA) 
follow-up program (Riess et al. 1999a), Riess et al. (1996a) show that 
the dispersion decreases from 0.52 mag to 0.12 mag after application of 
this "multi-color light-curve shape" (MLCS) method. The results from 
an expanded set of nearly 50 SNe la indicate that the dispersion de- 
creases from 0.44 mag to 0.15 mag (Figure 1). The resulting Hubble 
constant is 65 ± 2 (statistical) ±7 (systematic) km s _1 Mpc" 1 , with an 
additional systematic and zero-point uncertainty of ±5 km s _1 Mpc -1 . 
(Re-calibrations of the Cepheid distance scale, and other recent refine- 
ments, lead to a best estimate of Hq = 72 ± 8 km s _1 Mpc -1 , where the 
error bar includes both statistical and systematic uncertainties; Freed- 
man et al. 2001.) Saha et al. (2001) still argue for a longer distance 
scale, with Hq 60 km s _1 Mpc -1 , using different choices in measuring 
the supernova distances (Parodi et al. 2000) and the Cepheid distances; 
see the discussions by Gibson et al. (2000) and Jha (2002). 

Apart from the controversy regarding their true peak luminosity, low- 
redshift SNe la provide the best evidence that the Hubble flow is linear 
(Riess et al. 1996a), with deviations from linearity well-explained by 
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random peculiar velocities and the local flow field (Riess et al. 1997; 
Jha 2002). Phillips et al. (1999) and Riess et al. (1996a) further argue 
that the dust affecting SNe la is not of circumstellar origin, and show 
quantitatively that the extinction curve in external galaxies typically 
does not differ from that in the Milky Way (cf. Branch & Tammann 
1992, but see Tripp 1998). 
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Figure 1: Hubble diagrams for SNe la (A. G. Riess 2001, private com- 
munication) with velocities (km s _1 ) in the COBE rest frame on the 
Cepheid distance scale. The ordinate shows distance modulus, m — M 
(mag). Top: The objects are assumed to be standard candles and there 
is no correction for extinction; the result is a = 0.42 mag and Hq = 58±8 
km s _1 Mpc -1 . Bottom: The same objects, after correction for redden- 
ing and intrinsic differences in luminosity. The result is a = 0.15 mag 
and Hq = 65 ± 2 (statistical) km s _1 Mpc -1 , subject to changes in the 
zero-point of the Cepheid distance scale. (Indeed, the latest results with 
SNe la favor H = 72 km s" 1 Mpc" 1 .) 
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The advantage of systematically correcting the luminosities of SNe la 
at high redshifts rather than trying to isolate "normal" ones seems clear 
in view of evidence that the luminosity of SNe la may be a function of 
stellar population. If the most luminous SNe la occur in young stellar 
populations (e.g., Hamuy et al. 1996a, 2000; Branch, Romanishin, & 
Baron 1996; Ivanov, Hamuy, & Pinto 2000), then we might expect the 
mean peak luminosity of high-z SNe la to differ from that of a local 
sample. Alternatively, the use of Cepheids (Population I objects) to 
calibrate local SNe la can lead to a zero-point that is too luminous. 
On the other hand, as long as the physics of SNe la is essentially the 
same in young stellar populations locally and at high redshift, we should 
be able to adopt the luminosity correction methods (photometric and 
spectroscopic) found from detailed studies of low- z SNe la. 

In the past few years, many nearby SNe have been found by industri- 
ous amateur astronomers including R. Arbour, M. Armstrong, T. Boles, 
T. Puckett, M. Schwartz, and others. The "Nearby Supernova Factory" 
run by G. Aldering's team at the Lawrence Berkeley National Labo- 
ratory is also responsible for many discoveries, when it is conducting 
a search. However, about half of all reported nearby SNe in the last 
half-decade were discovered by my team's Lick Observatory Supernova 
Search (LOSS) conducted with the 0.76-m Katzman Automatic Imaging 
Telescope (KAIT; Li et al. 2000; Filippenko et al. 2001; Filippenko 2003; 
see http://astro.berkeley.edu/~bait/kait.html). During LOSS, CCD im- 
ages are taken of ^ 1000 galaxies per night and compared with KAIT 
"template images" obtained earlier; the templates are automatically sub- 
tracted from the new images and analyzed with computer software, and 
the SN candidates are flagged. The next day, undergraduate students at 
UC Berkeley examine all candidates, including weak ones, to eliminate 
star-like cosmic rays, asteroids, and other sources of false alarms. They 
also glance at some of the subtracted images to locate SNe that might be 
near bright, poorly subtracted stars or galactic nuclei. LOSS discovered 
20 SNe in 1998, 40 in 1999, 38 in 2000, 68 in 2001, 82 in 2002, and 95 
in 2003, making it by far the world's leading search for nearby SNe. 

Spectroscopic classifications of these and other low-redshift SNe are 
provided by a number of groups (see the I A U Circulars) , including our 
own observations with the Lick 3-m and Keck telescopes, the CfA SN 
monitoring campaign, the Texas/McDonald Observatory group, the Asi- 
ago team, the European Research and Training Network, the Australian 
National University team, etc. The most important objects were pho- 
tometrically monitored with KAIT through BVRI (and sometimes U) 
filters (e.g., Li et al. 2001a, 2003b; Modjaz et al. 2001; Leonard et al. 
2002a,b; Foley et al. 2003), and unfiltered follow-up observations (e.g., 
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Matheson et al. 2001) were made of almost all of them during the course 
of LOSS. Note that recently, KAIT has also been used to automatically 
monitor the optical afterglows of gamma-ray bursts (Li et al. 2003a, c; 
Matheson et al. 2003). 

This growing sample of well-observed SNe la allows us to more pre- 
cisely calibrate the distance determinations, as well as to look for correla- 
tions between the observed properties of the SNe and their environment 
(Hubble type of host galaxy, metallicity, stellar population, etc.). Jha 
(2002) present UBVRI follow-up observations from the CfA SN mon- 
itoring campaign of 44 nearby SNe la (including many discovered by 
KAIT), from which they derive "MLCS2k2," an updated light-curve 
fitting technique that includes U -band templates (critical for SN la ob- 
servations at high redshift). Recently, progress has also been made in 
the near-infrared, with Krisciunas et al. (2001, 2003, 2004) presenting 
JHK light curves of nearby SNe la. These have enabled Krisciunas, 
Phillips, & Suntzeff (2004) to construct near-infrared Hubble diagrams 
of SNe la, with significantly less effect from dust extinction along the line 
of sight. Their results suggest that SNe la are much closer to standard 
candles in the near-infrared than in the optical, with little dependence 
of light-curve shape on near-infrared luminosity. 

3. Cosmological Uses: High Redshifts 

These same techniques can be applied to construct a Hubble diagram 
with high-redshift SNe la, from which the value of qo = (Qm/2) — 
can be determined. With enough objects spanning a range of redshifts, 
we can measure Q,m and £l\ independently (e.g., Goobar & Perlmutter 
1995). Contours of peak apparent i?-band magnitude for SNe la at two 
redshifts have different slopes in the S7m - ^a plane, and the regions of 
intersection provide the answers we seek. 

3.1 The Search 

Based on the pioneering work of Norgaard-Nielsen et al. (1989), whose 
goal was to find SNe in moderate-redshift clusters of galaxies, the SCP 
(Perlmutter et al. 1995a, 1997) and the HZT (Schmidt et al. 1998) 
devised a strategy that almost guarantees the discovery of many faint, 
distant SNe la "on demand," during a predetermined set of nights. This 
"batch" approach to studying distant SNe allows follow-up spectroscopy 
and photometry to be scheduled in advance, resulting in a systematic 
study not possible with random discoveries. Most of the searched fields 
are equatorial, permitting follow-up from both hemispheres. The SCP 
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was the first group to convincingly demonstrate the ability to find SNe 
in batches. 

Our approach is simple in principle; see Schmidt et al. (1998) for de- 
tails, and for a description of our first high-redshift SN la (SN 1995K). 
Pairs of first-epoch images during the nights around new moon are ob- 
tained with wide-field cameras on large telescopes (e.g., the Big Through- 
put Camera on the CTIO 4-m Blanco telescope, and more recently the 
CTIO Mosaic II, the CFHT 8K and 12K mosaics, and Suprime Cam 
on Subaru), followed by second-epoch images 3-4 weeks later. (Pairs 
of images permit removal of cosmic rays, asteroids, and Kuiper-belt ob- 
jects.) These are compared immediately using well-tested software, and 
new SN candidates are identified in the second-epoch images (Figure 
2). Spectra are obtained as soon as possible after discovery to verify 
that the objects are SNe la and to determine their redshifts. Each team 
has already found over 200 SNe in concentrated batches, as reported in 
numerous IAU Circulars (e.g., Perlmutter et al. 1995b, 11 SNe with 
0.16 < z < 0.65; Suntzeff et al. 1996, 17 SNe with 0.09 < z < 0.84). 




Figure 2: Discovery image of SN 1997cj (28 April 1997), along with 
the template image and an HST image obtained subsequently. The net 
(subtracted) image is also shown. 

Intensive photometry of the SNe la commences within a few days af- 
ter procurement of the second-epoch images; it is continued throughout 
the ensuing and subsequent dark runs. In a few cases HST images are 
obtained. As expected, most of the discoveries are on the rise or near 
maximum brightness. When possible, the SNe are observed in filters that 
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closely match the redshifted B and V bands; this way, the K-corrections 
become only a second-order effect (Kim, Goobar, & Perhnutter 1996; 
Nugent, Kim, & Perlmutter 2002). Observations through two filters al- 
low us to apply reddening and luminosity corrections (Riess et al. 1996a; 
Hamuy et al. 1996a,b; Phillips et al. 1999). 

Although SNe in the magnitude range 22-22.5 can sometimes be spec- 
troscopically confirmed with 4-m class telescopes, the signal-to-noise ra- 
tios are low, even after several hours of integration. Certainly Keck, 
Gemini, the VLT, Subaru, or Magellan are required for the fainter ob- 
jects (22.5-24.5 mag). With the largest telescopes, not only can we 
rapidly confirm a substantial number of candidate SNe, but we can 
search for peculiarities in the spectra that might indicate evolution of 
SNe la with redshift. Moreover, high-quality spectra allow us to mea- 
sure the age of a SN: we have developed a method for automatically 
comparing the spectrum of a SN la with a library of spectra from many 
different epochs in the development of SNe la (Riess et al. 1997). Our 
technique also has great practical utility at the telescope: we can deter- 
mine the age of a SN "on the fly," within half an hour after obtaining its 
spectrum. This allows us to decide rapidly which SNe are best for subse- 
quent photometric follow-up, and we immediately alert our collaborators 
elsewhere. 

3.2 Results 

First, we note that the light curves of high-redshift SNe la are broader 
than those of nearby SNe la; the initial indications (Leibundgut et al. 
1996; Goldhaber et al. 1997), based on small numbers of SNe la, are 
amply confirmed with the larger samples (Goldhaber et al. 2001). Quan- 
titatively, the amount by which the light curves are "stretched" is con- 
sistent with a factor of 1 + z, as expected if redshifts are produced by the 
expansion of space rather than by "tired light" and other non-expansion 
hypotheses for the redshifts of objects at cosmological distances. [For 
non-standard cosmological interpretations of the SN la data, see Narlikar 
& Arp (1997) and Hoyle, Burbidge, & Narlikar (2000).] We also demon- 
strate this spectroscopically at the 2a confidence level for a single object: 
the spectrum of SN 1996bj (z = 0.57) evolved more slowly than those 
of nearby SNe la, by a factor consistent with 1 + z (Riess et al. 1997). 
Observations of SN 1997ex (z = 0.36) at three epochs conclusively verify 
the effects of time dilation: temporal changes in the spectra are slower 
than those of nearby SNe la by roughly the expected factor of 1.36. Al- 
though one might be able to argue that something other than universal 
expansion could be the cause of the apparent stretching of SN la light 
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curves at high redshifts, it is much more difficult to attribute apparently 
slower evolution of spectral details to an unknown effect. 

The formal value of CIm derived from SNe la has changed with time. 
The SCP published the initial result (Perlmutter et al. 1995a), based on 
a single object, SN 1992bi at z = 0.458: U M = 0.2 ± 0.6 ± 1.1 (assuming 
that = 0). The SCP's analysis of their first seven objects (Perlmutter 
et al. 1997) suggested a much larger value of = 0.88±0.6 (if S7a = 0) 
or Qm = 0.94 ± 0.3 (if J^totai = 1)- Such a high-density universe seemed 
at odds with other, independent measurements of VLm at that time. 
However, with the subsequent inclusion of just one more object, SN 
1997ap at z = 0.83 (the highest then known for a SN la; Perlmutter et 
al. 1998), their estimates were revised back down to Q,m = 0.2 ± 0.4 if 
f&A = 0, and CIm = 0.6 =L 0.2 if O to tai = 1> the apparent brightness of 
SN 1997ap had been precisely measured with HST, so it substantially 
affected the best fits. 

Meanwhile, the HZT published (Garnavich et al. 1998a) an analysis 
of four objects (three of them observed with HST), including SN 1997ck 
at z = 0.97, at that time a redshift record, although they cannot be 
absolutely certain that the object was a SN la because the spectrum 
is too poor. From these data, the HZT derived that Hm = —0.1 ± 
0.5 (assuming £l\ = 0) and Qm = 0.35 ± 0.3 (assuming Sltotai = 1), 
inconsistent with the high initially found by Perlmutter et al. (1997) 
but consistent with the revised estimate in Perlmutter et al. (1998). An 
independent analysis of 10 SNe la using the "snapshot" distance method 
(with which conclusions are drawn from sparsely observed SNe la) gave 
quantitatively similar conclusions (Riess et al. 1998a). However, none 
of these early data sets carried the statistical discriminating power to 
detect cosmic acceleration. 

The SCP's next results were announced at the 1998 January AAS 
meeting in Washington, DC. A press conference was scheduled, with the 
stated purpose of presenting and discussing the then-current evidence 
for a low-^Af universe as published by Perlmutter et al. (1998; SCP) 
and Garnavich et al. (1998a; HZT). When showing the SCP's Hubble 
diagram for SNe la, however, Perlmutter also pointed out tentative ev- 
idence for acceleration). He said that the conclusion was uncertain, and 
that the data were equally consistent with no acceleration; the system- 
atic errors had not yet been adequately assessed. Essentially the same 
conclusion was given by the SCP in February 1998 during their talks at 
a major conference on dark matter held near Los Angeles (Goldhaber & 
Perlmutter 1998). 

Although it chose not to reveal them at the same 1998 January AAS 
meeting, the HZT already had similar, tentative evidence for acceleration 
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in their own SN la data set. The HZT continued to perform numerous 
checks of their data analysis and interpretation, including fairly thorough 
consideration of various possible systematic effects. Unable to find any 
significant problems, even with the possible systematic effects, the HZT 
reported detection of a nonzero value for f^A (based on 16 high-z SNe la) 
at the Los Angeles dark matter conference in February 1998 (Filippenko 
& Riess 1998), and soon thereafter submitted a formal paper that was 
published in September 1998 (Riess et al. 1998b). 




Figure 3 (left): The upper panel shows the Hubble diagram for the low- 
z and high-z HZT SN la sample with MLCS distances; see Riess et 
al. (1998b). Overplotted are three world models: "low" and "high" Qm 
with = 0, and the best fit for a flat universe (Qm = 0.28, f^A = 0.72). 
The bottom panel shows the difference between data and models from 
the = 0.20, S7a = prediction. Only the 10 best-observed high-z 
SNe la are shown. The average difference between the data and the 
VLm = 0.20, Q\ = prediction is ~ 0.25 mag. 

Figure 4 (right): The HZT's combined constraints from SNe la (left) 
and the position of the first acoustic peak of the cosmic microwave 
background (CMB) angular power spectrum, based on data available in 
mid-1998; see Garnavich et al. (1998b). The contours mark the 68.3%, 
95.4%, and 99.7% enclosed probability regions. Solid curves correspond 
to results from the MLCS method; dotted ones are from the Ami5(i?) 
method; all 16 SNe la in Riess et al. (1998b) were used. 
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The HZT's original Hubble diagram for the 10 best-observed high-z 
SNe la is given in Figure 3. With the MLCS method applied to the 
full set of 16 SNe la, the HZT's formal results were Q, M = 0.24 ± 0.10 
if fitotai = 1) or = —0.35 ± 0.18 (unphysical) if f^A = 0. If one 
demanded that = 0.2, then the best value for Oa was 0.66 ± 0.21. 
These conclusions did not change significantly when only the 10 best- 
observed SNe la were used (Figure 3; Qm = 0.28 ± 0.10 if fitotai = !)■ 

Another important constraint on the cosmological parameters could 
be obtained from measurements of the angular scale of the first acoustic 
peak of the CMB (e.g., Zaldarriaga, Spergel, & Seljak 1997; Eisenstein, 
Hu, & Tegmark 1998); the SN la and CMB techniques provide nearly 
complementary information. A stunning result was already available 
by mid-1998 from existing measurements (e.g., Hancock et al. 1998; 
Lineweaver &; Barbosa 1998): the HZT's analysis of the SN la data in 
Riess et al. (1998b) demonstrated that Q M + = 0.94 ± 0.26 (Figure 
4), when the SN and CMB constraints were combined (Garnavich et al. 
1998b; see also Lineweaver 1998, Efstathiou et al. 1999, and others). 

Somewhat later (June 1999), the SCP published almost identical re- 
sults (Perlmutter et al. 1999), implying an accelerating expansion of the 
Universe, based on an essentially independent set of 42 high-z SNe (al- 
though a substantial number of these were later discarded because they 
are not clearly SNe la; Knop et al. 2003). The Perlmutter et al. (1999) 
data, together with those of the HZT, are shown in Figure 5, and the 
corresponding confidence contours in the f^A vs. plane are given in 
Figure 6. This incredible agreement suggested that neither group had 
made a large, simple blunder; if the result was wrong, the reason must 
be subtle. Had there been only one team working in this area, it is likely 
that far fewer astronomers and physicists throughout the world would 
have taken the result seriously. 

Moreover, already in 1998-1999 there was tentative evidence that the 
"dark energy" driving the accelerated expansion was indeed consistent 
with the cosmological constant, A. If A dominates, then the equation-of- 
state parameter of the dark energy should bew = —1, where the pressure 
(P) and density (p) are related according to w = P/(pc 2 ). Garnavich et 
al. (1998b) and Perlmutter et al. (1999) were able to set an interesting 
limit, w ^ — 0.60 at the 95% confidence level. However, more high- 
quality data at z « 0.5 are needed to narrow the allowed range, in order 
to test other proposed candidates for dark energy such as various forms 
of "quintessence" (e.g., Caldwell, Dave, & Steinhardt 1998). 

Although the CMB results appeared reasonably persuasive in 1998- 
1999, one could argue that fluctuations on different scales had been 
measured with different instruments, and that subtle systematic ef- 
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fects might lead to erroneous conclusions. These fears were dispelled 
only 1-2 years later, when the more accurate and precise results of the 
BOOMERANG collaboration were announced (de Bernardis et al. 2000, 
2002). Shortly thereafter the MAXIMA collaboration distributed their 
very similar findings (Hanany et al. 2000; Balbi et al. 2000; Netterfield 
et al. 2002; see also the TOCO, DASI, and many other measurements). 
Figure 6 illustrates that the CMB measurements tightly constrain fitotai 
to be close to unity; we appear to live in a flat universe, in agreement 
with most inflationary models for the early Universe! Combined with 
the SN la results, the evidence for nonzero f^A was fairly strong. Making 
the argument even more compelling was the fact that various studies of 
clusters of galaxies (see summary by Bahcall et al. 1999) showed that 
m 0.3, consistent with the results in Figures 4 and 6. Thus, a "con- 
cordance cosmology" had emerged: £Im ~ 0.3, f^A ~ 0.7 — consistent 
with what had been suspected some years earlier by Ostriker & Stein- 
hardt (1995; see also Carroll, Press, k Turner 1992). 




0.0 0.5 1.0 1.5 2.0 2.5 

Figure 5 (left): As in Figure 3, but this time including both the HZT 
(Riess et al. 1998b) and SCP (Perlmutter et al. 1999) samples of low- 
redshift and high-redshift SNe la. Overplotted are three world models: 
£Im = 0.3 and 1.0 with f^A = 0, and a flat universe (^totai = 1-0) with 
f^A = 0.7. The bottom panel shows the difference between data and 
models from the Qm = 0.3, f^A = prediction. 

Figure 6 (right): The combined constraints from SNe la (see Figure 5) 
and the position of the first acoustic peak of the CMB angular power 
spectrum, based on BOOMERANG and MAXIMA data. The contours 
mark the 68.3%, 95.4%, and 99.7% enclosed probability regions deter- 
mined from the SNe la. According to the CMB, Ototai ~ 1-0. 
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Yet another piece of evidence for a nonzero value of A was provided by 
the Two-Degree Field Galaxy Redshift Survey (2dFGRS; Peacock et al. 
2001; Percival et al. 2001; Efstathiou et al. 2002). Combined with the 
CMB maps, their results are inconsistent with a universe dominated by 
gravitating dark matter. Again, the implication is that about 70% of the 
mass-energy density of the Universe consists of some sort of dark energy 
whose gravitational effect is repulsive. In 2003, results from the first 
year of Wilkinson Microwave Anisotropy Probe (WMAP) observations 
appeared; together with the 2dFGRS constraints, they confirmed and 
refined the concordance cosmology (£Im = 0.27, ^a = 0.73, ^baryon = 
0.044, H = 71 ± 4 km s" 1 Mpc -1 ; Spergel et al. 2003, and references 
therein). Recent evidence for dark energy and an accelerating universe 
has also come from 2-3cr detections of the integrated Sachs- Wolfe effect 
(Afshordi, Loh, & Strauss 2004; Boughn & Crittenden 2004; Fosalba et 
al. 2003; Nolta et al. 2004; Scranton et al. 2004). 

Meanwhile, the SN la measurements were becoming more numerous 
and of higher quality. For the HZT, the new (Fall 1999) sample of high- 
redshift SNe la presented by Tonry et al. (2003), analyzed with methods 
distinct from (but similar to) those used previously, confirmed the result 
of Riess et al. (1998b) and Perlmutter et al. (1999) that the expansion of 
the Universe is accelerating. By combining all of the data sets available 
at that time, Tonry et al. (2003) were able to use about 200 SNe la, 
obtaining an incredibly firm detection of f^A > 0. They placed the 
following constraints on cosmological quantities: (1) If the equation-of- 
state parameter of the dark energy is w = —1, then i^o^o = 0.96 ± 0.04, 
and Q\ — 1AQ.M = 0.35 ± 0.14. (2) Including the constraint of a flat 
universe, they find that Qm = 0.28±0.05, independent of any large-scale 
structure measurements. (3) Adopting a prior based on the 2dFGRS 
constraint on £Im (Percival et al. 2001) and assuming a flat universe, 
they derive that -1.48 < w < -0.72 at 95% confidence. (4) Adopting 
the 2dFGRS results, they find Vt M = 0.28 and n A = 0.72, independent of 
any assumptions about ^totai- These constraints are similar in precision 
and in value to conclusions reported using Year 1 of WMAP (Spergel et 
al. 2003), also in combination with the 2dFGRS. 

Even more recently, Barris et al. (2004) describe results from the 2001 
HZT campaign. 23 SNe with z = 0.34-1.03 were studied, 9 of which are 
unambiguously SNe la. The sample includes 15 objects at z > 0.7, 
doubling the number of published SNe at these redshifts. Under the 
assumption that fitotai = 1> they determine best-fit values of (Q^, Q\) 
= (0.33, 0.67), using 22 SNe from the 2001 HZT survey together with 
172 SNe from Tonry et al. (2003) satisfying z > 0.01 and Ay < 0.5 mag. 
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The SCP's high-quality HST data set on 11 SNe la in the redshift 
range 0.36-0.86, recently published by Knop et al. (2003), indepen- 
dently confirms the apparent acceleration of the Universe. They were 
able to measure accurate colors of the SNe, providing better host-galaxy 
extinction estimates than had been possible in the past for individual 
objects. Thus, there was no need to make any assumptions or priors 
on the parent distribution of extinction values, E(B — V), unlike the 
case in the initial analyses (Riess et al. 1998; Perlmutter et al. 1999; 
but see Sullivan et al. 2003, discussed below). Their extinction mea- 
surements do not show evidence for anomalously blue SNe la at high 
redshifts (or for a preponderance of negative E(B — V) values for high- 
redshift SNe la), contrary to some earlier suspicions (Falco et al. 1999; 
Leibundgut 2001). Knop et al. (2003) find that dark energy is required 
with a probability exceeding 99%, consistent with previous studies that 
had made assumptions about the distribution of extinctions or that had 
used low-extinction subsets of SNe la. In a flat universe with a con- 
stant dark energy equation-of-state parameter w = — 1, they find that 
Q A = 0.75 (+0.06, -0.07) ± 0.04, where the first quoted uncertainties 
are statistical and the second are identified systematics. Moreover, in 
a flat universe with w constant in time, their SNe la data show that 
w = -1.05(+0.15, -0.20) ± 0.09, consistent with A (rather than with 
some quintessence models; see below). 

The dynamical age of the Universe can be calculated from the cos- 
mological parameters. In an empty Universe with no cosmological con- 
stant, the dynamical age is simply the inverse of the Hubble constant, 
to = -f^o" 1 ; there is no deceleration. In the late-1990s, SNe la gave 
H = 65 ± 7 km s _1 Mpc" 1 , and a Hubble time of 15.1 ±1.6 Gyr. For a 
more complex cosmology, integrating the velocity of the expansion from 
the current epoch (z = 0) to the beginning (z = oo) yields an expression 
for the dynamical age. As shown in detail by Riess et al. (1998b), by 
mid-1998 the HZT had obtained a value of 14.2+J;^ Gyr (with H = 65) 
using the likely range for (Om,^a) that they measured. (The preci- 
sion was so high because their experiment was sensitive to roughly the 
difference between and f^A, and the dynamical age also varies in 
approximately this way.) Including the systematic uncertainty of the 
Cepheid distance scale, which may be up to 10%, a reasonable estimate 
of the dynamical age was 14.2 ±1.7 Gyr (Riess et al. 1998b). Again, the 
SCP's result was very similar (Perlmutter et al. 1999), since it was based 
on nearly the same derived values for the cosmological parameters. The 
most recent results, reported by Tonry et al. (2003; see also Knop et al. 
2003) and adopting H = 72 ±8 km s" 1 Mpc -1 , give a dynamical age of 
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13.1 ± 1.5 Gyr for the Universe — again, in agreement with the WMAP 
result of 13.7 ± 0.2 Gyr. 

This expansion age is also consistent with ages determined from var- 
ious other techniques such as the cooling of white dwarfs (Galactic disk 
> 9.5 Gyr; Oswalt et al. 1996), radioactive dating of stars via the tho- 
rium and europium abundances (15.2 ±3.7 Gyr; Cowan et al. 1997), and 
studies of globular clusters (10-15 Gyr, depending on whether Hipparcos 
parallaxes of Cepheids are adopted; Gratton et al. 1997; Chaboyer et 
al. 1998). The ages of the oldest stars no longer seem to exceed the ex- 
pansion age of the Universe; the long-standing "age crisis" has evidently 
been resolved. 

4. Discussion 

Although the convergence of different methods on the same answer is 
reassuring, and suggests that the concordance cosmology is correct, it is 
important to vigorously test each method to make sure it is not leading 
us astray. Moreover, only through such detailed studies will the accuracy 
and precision of the methods improve, allowing us to eventually set 
better constraints on the equation-of-state parameter, w. Here I discuss 
the systematic effects that could adversely affect the SN la results. 

High-redshift SNe la are observed to be dimmer than expected in 
an empty Universe (i.e., SIm = 0) with no cosmological constant. At 
z ~ 0.5, where the SN la observations have their greatest leverage on 
A, the difference between an = 0.3 (Oa = 0) universe and a flat 
universe with £l\ = 0.7 is only about 0.25 mag. Thus, we need to find out 
if chemical abundances, stellar populations, selection bias, gravitational 
lensing, or grey dust can have an effect this large. Although both the 
HZT and SCP had considered many of these potential systematic effects 
in their original discovery papers (Riess et al. 1998b; Perlmutter et al. 
1999), and had shown with reasonable confidence that obvious ones were 
not greatly affecting their conclusions, it was of course possible that they 
were wrong, and that the data were being misinterpreted. 

4.1 Gravitational Lensing 

The magnification and demagnification of light by large-scale struc- 
ture can alter the observed magnitudes of high-redshift SNe (e.g., Kan- 
towski, Vaughan, & Branch 1995; Kantowski 1998). The effect of weak 
gravitational lensing on our analysis has been quantified by Wambsganss, 
Cen, & Ostriker (1998) and summarized by Schmidt et al. (1998). SN la 
light will, on average, be demagnified by 0.5% at z = 0.5 and by 1% 
at z = 1 in a Universe with a non- negligible cosmological constant. Al- 
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though the sign of the effect is the same as the influence of a cosmological 
constant, the size of the effect is small enough to be ignored. 

Holz (1998), Holz & Wald (1998), and Kantowski (1998) have cal- 
culated the weak lensing effects on supernova light from ordinary mat- 
ter which is not smoothly distributed in galaxies but rather clumped 
into stars (i.e., dark matter contained in massive compact halo objects). 
With this scenario, microlensing becomes a more important effect, fur- 
ther decreasing the typical observed supernova luminosities at z = 0.5 
by 0.02 mag for VLm = 0.2. Even if most ordinary matter were contained 
in compact objects, this effect would not be large enough to reconcile 
the SN la distances with the influence of ordinary matter alone in a 
decelerating universe. Barris et al. (2004), for example, show that use 
of the "empty-beam model" for gravitational lensing does not eliminate 
the need for > 0. However, gravitational lensing will certainly need 
to be taken into account when making a precise measurement of w. 

With a very large sample (200-1000) of high-redshift SNe la, and 
accurate photometry, it should be possible to quantify the effects of 
clumped matter. Light from some SNe la should be strongly amplified 
by the presence of intervening matter, while the vast majority will be 
deamp lined (e.g., Barber et al. 2000). The distribution of amplification 
factors can be used to determine the type of dark matter most prevalent 
in the Universe (compact objects, or smoothly distributed). 

4.2 Evolution 

Perhaps the most obvious possible systematic effect is evolution of 
SNe la over cosmic time, due to changes in metallicity, progenitor mass, 
or some other factor. If the peak luminosity of SNe la were lower at 
high redshift, then the case for S7a > would weaken. Conversely, if 
the distant explosions are more powerful, then the case for acceleration 
strengthens. Theorists are not yet sure what the sign of the effect will be, 
if it is present at all; different assumptions lead to different conclusions 
(Hoflich, Wheeler, & Thielemann 1998; Umeda et al. 1999; Nomoto et 
al. 2000; Yungelson & Livio 2000). 

Of course, it is extremely difficult, if not effectively impossible, to 
obtain an accurate, independent measure of the peak luminosity of high- 
z SNe la, and hence to directly test for luminosity evolution. However, 
we can more easily determine whether other observable properties of low- 
z and high-z SNe la differ. If they are all the same, it is more probable 
that the peak luminosity is constant as well — but if they differ, then 
the peak luminosity might also be affected (e.g., Hoflich et al. 1998). 
Drell, Loredo, & Wasserman (2000), for example, argue that there are 
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reasons to suspect evolution, because the average properties of existing 
samples of high- 2; and low-z SNe la seem to differ (e.g., the high- 2; SNe la 
are more uniform). 

The local sample of SNe la displays a weak correlation between light- 
curve shape and host-galaxy type, in the sense that the most luminous 
SNe la with the broadest light curves only occur in late- type galaxies. 
Both early-type and late-type galaxies provide hosts for dimmer SNe la 
with narrower light curves (Hamuy et al. 1996a). The mean luminosity 
difference for SNe la in late-type and early-type galaxies is ~ 0.3 mag. 
In addition, the SN la rate per unit luminosity is almost twice as high 
in late-type galaxies as in early-type galaxies at the present epoch (Cap- 
pellaro et al. 1997). These results may indicate an evolution of SNe la 
with progenitor age. Possibly relevant physical parameters are the mass, 
metallicity, and C/O ratio of the progenitor (Hdflich et al. 1998). 

We expect that the relation between light-curve shape and peak lu- 
minosity that applies to the range of stellar populations and progenitor 
ages encountered in the late-type and early-type hosts in our nearby 
sample should also be applicable to the range encountered in our distant 
sample. In fact, the age range for SN la progenitors in the nearby sample 
is likely larger than the change in mean progenitor age over the 4-6 Gyr 
lookback time to the high-z sample. Thus, to first order at least, our 
local sample should correct the distances for progenitor or age effects. 

We can place empirical constraints on the effect that a change in 
the progenitor age would have on our SN la distances by comparing 
subsamples of low-redshift SNe la believed to arise from old and young 
progenitors. In the nearby sample, the mean difference between the 
distances for the early-type hosts (8 SNe la) and late- type hosts (19 
SNe la), at a given redshift, is 0.04 ± 0.07 mag from the MLCS method. 
This difference is consistent with zero. Even if the SN la progenitors 
evolved from one population at low redshift to the other at high redshift, 
we still would not explain the surplus in mean distance of 0.25 mag over 
the = prediction. 

Moreover, in a major study of high-redshift SNe la as a function of 
galaxy morphology, the SCP found no clear differences (except for the 
amount of scatter; see §5.2) between the cosmological results obtained 
with SNe la in late- type and early- type galaxies (Sullivan et al. 2003). 
Similarly, a study of the host galaxies of high-redshift SNe la (0.42 < z < 
1.06) done by the HZT (Williams et al. 2003) found no clear evidence for 
correlations between host-galaxy properties and the residuals of distance 
measurements from cosmological fits. Some of the correlations between 
SN la properties and host-galaxy type seen at low redshift also appear 
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to be present at high redshift, again supporting the current practice of 
extrapolating properties of the nearby SN la population to high redshifts. 

It is also reassuring that initial comparisons suggest that high-z SN la 
spectra appear similar to those observed at low redshift. For example, 
the spectral characteristics of SN 1998ai (z = 0.49) appear to be essen- 
tially indistinguishable from those of normal low-z SNe la; see Figure 
7. In fact, the most obviously discrepant spectrum in this figure is the 
second one, that of SN 1994B (z = 0.09); it is intentionally included as a 
"decoy" that illustrates the degree to which even the spectra of nearby, 
relatively normal SNe la can vary. Nevertheless, it is important to note 
that a dispersion in luminosity (perhaps 0.2 mag) exists even among the 
other, more normal SNe la shown in Figure 7; thus, our spectra of SN 
1998ai and other high-z SNe la are not yet sufficiently good for indepen- 
dent, precise determinations of peak luminosity from spectral features 
(Nugent et al. 1995). Many of them, however, are sufficient for ruling 
out other SN types (Figure 8), or for identifying gross peculiarities such 
as those shown by SNe 1991T and 1991bg; see Coil et al. (2000). 

We can help verify that the SNe at z ss 0.5 being used for cosmol- 
ogy do not belong to a subluminous population of SNe la by examining 
rest-frame /-band light curves. Normal, nearby SNe la show a pro- 
nounced second maximum in the / band about a month after the first 
maximum and typically about 0.5 mag fainter (e.g., Ford et al. 1993; 
Suntzeff 1996). Subluminous SNe la, in contrast, do not show this sec- 
ond maximum, but rather follow a linear decline or show a muted second 
maximum (Filippenko et al. 1992a). As discussed by Riess et al. (2000), 
tentative evidence for the second maximum is seen from the HZT's ex- 
isting J-band (rest-frame /-band) data on SN 1999Q (z = 0.46); see 
Figure 10. Additional tests with spectra and near-infrared light curves 
are currently being conducted. 

Another way of using light curves to test for possible evolution of 
SNe la is to see whether the rise time (from explosion to maximum 
brightness) is the same for high-redshift and low-redshift SNe la; a differ- 
ence might indicate that the peak luminosities are also different (Hoflich 
et al. 1998). Riess et al. (1999c) measured the risetime of nearby SNe la, 
using data from KAIT, the Beijing Astronomical Observatory (BAO) SN 
search, and a few amateur astronomers. Though the exact value of the 
risetime is a function of peak luminosity, for typical low-redshift SNe la 
it is 20.0 ±0.2 days. Riess et al. (1999b) pointed out that this differs by 
5.8cr from the preliminary risetime of 17.5 ± 0.4 days reported in confer- 
ences by the SCP (Goldhaber et al. 1998a,b; Groom 1998). However, 
more thorough analyses of the SCP data (Aldering, Knop, & Nugent 
2000; Goldhaber et al. 2001) show that the high-redshift uncertainty of 
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±0.4 days that the SCP originally reported was much too small because 
it did not account for systematic effects. The revised discrepancy with 
the low-redshift risetime is about 2a or less. Thus, the apparent differ- 
ence in risetimes might be insignificant. Even if the difference is real, 
however, its relevance to the peak luminosity is unclear; the light curves 
may differ only in the first few days after the explosion, and this could 
be caused by small variations in conditions near the outer part of the 
exploding white dwarf that are inconsequential at the peak. 
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Figure 7 (left): Spectral comparison (in fx) of SN 1998ai (z = 0.49; Keck 
spectrum) with low-redshift (z < 0.1) SNe la at a similar age (~ 5 days 
before maximum brightness), from Riess et al. (1998b). The spectra 
of the low-redshift SNe la were resampled and convolved with Gaussian 
noise to match the quality of the spectrum of SN 1998ai. Overall, the 
agreement in the spectra is excellent, tentatively suggesting that distant 
SNe la are physically similar to nearby SNe la. SN 1994B (z = 0.09) 
differs the most from the others, and was included as a "decoy." 
Figure 8 (right): Heavily smoothed spectra of two high-z SNe (SN 1999ff 
at z = 0.455 and SN 1999fv at z = 1.19; quite noisy below ~3500 A) are 
presented along with several low-z SN la spectra (SNe 1989B, 1992A, 
and 1981B), a SN lb spectrum (SN 1993 J), and a SN Ic spectrum (SN 
19941); see Filippenko (1997b) for a discussion of spectra of various types 
of SNe. The date of the spectra relative to 5-band maximum is shown 
in parentheses after each object's name. Specific features seen in SN 
1999ff and labeled with a letter are discussed by Coil et al. (2000). This 
comparison shows that the two high-z SNe are most likely SNe la. 
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Although there are no clear signs that cosmic evolution of SNe la 
seriously compromises our results, it is wise to remain vigilant for possi- 
ble problems. At low redshifts, for example, we already know that some 
SNe la don't conform with the correlation between light curve shape and 
luminosity. SN 2000cx in the SO galaxy NGC 524, for example, has light 
curves that cannot be fit well by any of the fitting techniques currently 
available (Li et al. 2001a; Filippenko 2003); see Figure 9. Its late-time 
color is remarkably blue, inconsistent with the homogeneity described 
by Phillips et al. (1999). The spectral evolution of SN 2000cx is pe- 
culiar as well (Li et al. 2001a; Branch et al. 2004a): the photosphere 
appears to have remained hot for a long time, and both iron-peak and 
intermediate-mass elements move at very high velocities. 
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Figure 9: The MLCS fit (Riess et al. 1998b; left panel) and the stretch 
method fit (Perlmutter et al. 1999; right panel) for SN 2000cx. The 
MLCS fit is the worst we had ever seen through year 2000. For the 
stretch method fit, the solid line is the fit to all the data points from 
t = —8 to 32 days, the dash-dotted line uses only the premaximum 
datapoints, and the dashed line only the postmaximum datapoints. The 
three fits give very different stretch factors. From Li et al. (2001a). 



An even more peculiar object is SN 2002cx (Li et al. 2003; Fil- 
ippenko 2003; Branch et al. 2004b). It is spectroscopically bizarre, 
with extremely low expansion velocities and almost no evidence for 
intermediate-mass elements. The nebular phase was reached incredibly 
soon after maximum brightness, despite the low velocity of the ejecta, 
suggesting that the ejected mass is small. SN 2002cx was subluminous 
by ~ 2 mag at all optical wavelengths relative to normal SNe la, despite 
the early-time spectroscopic resemblance to the somewhat overluminous 
SN 1991T. The i?-band and /-band light curves of SN 2002cx are com- 
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pletely unlike those of normal SNe la. No existing theoretical model 
successfully explains all observed aspects of SN 2002cx, though 3D de- 
flagration models may be best. If there are more strange beasts like SNe 
2000cx and 2002cx at high redshifts than at low redshifts, systematic 
errors may creep into the analysis of high-z SNe la. 

4.3 Extinction 

Our SN la distances have the important advantage of including cor- 
rections for interstellar extinction occurring in the host galaxy and the 
Milky Way. Extinction corrections based on the relation between SN la 
colors and luminosity improve distance precision for a sample of nearby 
SNe la that includes objects with substantial extinction (Riess et al. 
1996a; Phillips et al. 1999); the scatter in the Hubble diagram is much 
reduced. Moreover, the consistency of the measured Hubble flow from 
SNe la with late-type and early-type hosts (see §5.1) shows that the 
extinction corrections applied to dusty SNe la at low redshift do not al- 
ter the expansion rate from its value measured from SNe la in low-dust 
environments. 

In practice, the high-redshift SNe la generally appear to suffer very 
little extinction; their B — V colors at maximum brightness are normal, 
suggesting little color excess due to reddening. The most detailed avail- 
able study is that of the SCP (Sullivan et al. 2003): they found that the 
scatter in the Hubble diagram is minimal for SNe la in early-type host 
galaxies, but increases for SNe la in late-type galaxies. Moreover, on 
average the SNe in late- type galaxies are slightly fainter (by 0.14 ± 0.09 
mag) than those in early-type galaxies. Finally, at peak brightness the 
colors of SNe la in late-type galaxies are marginally redder than those 
in early-type galaxies. Sullivan et al. (2003) conclude that extinction by 
dust in the host galaxies of SNe la is one of the major sources of scat- 
ter in the high-redshift Hubble diagram. By restricting their sample to 
SNe la in early- type host galaxies (presumably with minimal extinction), 
they obtain a very tight Hubble diagram that suggests a nonzero value 
for at the 5a confidence level, under the assumption that fitotai = 1- 
In the absence of this assumption, SNe la in early-type hosts still imply 
that S7a > at nearly the 98% confidence level. The results for JIa with 
SNe la in late-type galaxies are quantitatively similar, but statistically 
less secure because of the larger scatter. 

Riess, Press, & Kirshner (1996b; see also Phillips et al. 1999) found 
indications that the Galactic ratios between selective absorption and 
color excess are similar for host galaxies in the nearby (z < 0.1) Hubble 
flow. Yet, what if these ratios changed with lookback time (e.g., Aguirre 
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1999a)? Could an evolution in dust-grain size descending from ancestral 
interstellar "pebbles" at higher redshifts cause us to underestimate the 
extinction? Large grains would not imprint the reddening signature of 
typical interstellar extinction upon which our corrections necessarily rely. 
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Figure 10 (left): Rest-frame /-band (observed J-band) light curve of 
SN 1999Q (z = 0.46, 5 solid points; Riess et al. 2000), and the /-band 
light curves of several nearby SNe la. Subluminous SNe la exhibit a less 
prominent second maximum than do normal SNe la. 

Figure 11 (right): Color (B - I) and color excess (Eb-i) for SN 1999Q 
and different dust models (Riess et al. 2000). The data are most con- 
sistent with no dust and f^A > 0. 

However, viewing our SNe through such gray interstellar grains would 
also induce a dispersion in the derived distances. Using the results of 
Hatano, Branch, & Deaton (1998), Riess et al. (1998b) estimate that 
the expected dispersion would be 0.40 mag if the mean gray extinction 
were 0.25 mag (the value required to explain the measured MLCS dis- 
tances without a cosmological constant). This is significantly larger than 
the 0.21 mag dispersion observed in the high-redshift MLCS distances. 
Furthermore, most of the observed scatter is already consistent with the 
estimated statistical errors, leaving little to be caused by gray extinction. 
Nevertheless, if we assumed that all of the observed scatter were due to 
gray extinction, the mean shift in the SN la distances would be only 
0.05 mag. With the existing observations, it is difficult to rule out this 
modest amount of gray interstellar extinction. 

Gray intergalactic extinction could dim the SNe without either tell- 
tale reddening or dispersion, if all lines of sight to a given redshift had 
a similar column density of absorbing material. The component of the 
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intergalactic medium with such uniform coverage corresponds to the gas 
clouds producing Lyman-a forest absorption at low redshifts. These 
clouds have individual H I column densities less than about 10 15 cm" 2 
(Bahcall et al. 1996). However, they display low metallicities, typically 
less than 10% of solar. Gray extinction would require larger dust grains 
which would need a larger mass in heavy elements than typical interstel- 
lar grain size distributions to achieve a given extinction. It is possible 
that large dust grains are blown out of galaxies by radiation pressure, 
and are therefore not associated with Lyman-a clouds (Aguirre 1999b). 

But even the dust postulated by Aguirre ( 1999a, b) and Aguirre & 
Haiman (1999) is not completely gray, having a size of about 0.1 ^m. We 
can test for such nearly gray dust by observing high-redshift SNe la over 
a wide wavelength range to measure the color excess it would introduce. 
UA V = 0.25 mag, then E(U-I) and E(B-I) should be 0.12-0.16 mag 
(Aguirre 1999a,b). If, on the other hand, the 0.25 mag faintness is due 
to A, then no such reddening should be seen. This effect is measurable 
using proven techniques; so far, with just one SN la (SN 1999Q; Figure 
11), our results favor the no-dust hypothesis to better than 2a (Riess et 
al. 2000). More work along these lines is in progress. 

4.4 Early Deceleration of the Universe 

Suppose, however, that for some reason the dust is very gray, or our 
color measurements are not sufficiently precise to rule out Aguirre's (or 
other) dust. Or, perhaps some other astrophysical systematic effect is 
fooling us, such as possible evolution of the white dwarf progenitors 
(e.g., Hoflich et al. 1998; Umeda et al. 1999), or gravitational lensing 
(Wambsganss et al. 1998). The most decisive test to distinguish between 
A and cumulative systematic effects is to examine the deviation of the 
observed peak magnitude of SNe la from the magnitude expected in the 
1ow-$7m ; zero- A model. If A is positive, the deviation should actually 
begin to decrease at z « 1; we will be looking so far back in time that 
the A effect becomes small compared with £Im, and the Universe is 
decelerating at that epoch. If, on the other hand, a systematic bias such 
as gray dust or evolution of the white dwarf progenitors is the culprit, 
we expect that the deviation of the apparent magnitude will continue 
growing, unless the systematic bias is set up in such an unlikely way as 
to mimic the effects of A (Drell et al. 2000). A turnover, or decrease of 
the deviation of apparent magnitude at high redshift, can be considered 
almost like a "smoking gun" of A (or, more generally, of a dark-energy 
component whose density does not change much with time). 
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In a wonderful demonstration of good luck and hard work, Riess et al. 
(2001) report on HST observations of a probable SN la at z ~ 1.7 (SN 
1997ff, the most distant SN ever observed) that suggest the expected 
turnover is indeed present, providing a tantalizing glimpse of the epoch 
of deceleration. (See also Bemtez et al. 2002, which corrects the observed 
magnitude of SN 1997ff for gravitational lensing by foreground galaxies.) 
SN 1997ff was discovered by Gilliland & Phillips (1998) in a repeat HST 
observation of the Hubble Deep Field-North, and it was serendipitously 
monitored in the infrared with ffi)T/NICMOS. The peak apparent SN 
brightness is consistent with that expected in the decelerating phase of 
the concordance cosmological model, Om ~ 0.3, rs 0.7 (Figure 12). 
It is inconsistent with gray dust or simple luminosity evolution, when 
combined with the data for SNe la at z ~ 0.5. 
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Figure 12: Hubble diagram for SNe la relative to an empty universe 
(f2 = 0) compared with cosmological and astrophysical models (Riess et 
al. 2001). Low-redshift SNe la are from Hamuy et al. (1996a) and Riess 
et al. (1999a). The magnitude of SN 1997ff at z = 1.7 has been corrected 
for gravitational lensing (Bemtez et al. 2002). The measurements of 
SN 1997ff are inconsistent with astrophysical effects that could mimic 
previous evidence for an accelerating universe from SNe la at z ~ 0.5. 



On the other hand, it was wise to remain cautious at the time: the 
error bars are large, and it is always possible that we are being fooled by 
this single object. The HZT and SCP thus started programs to find and 
measure more SNe la at such high redshifts {z > 1). One promising- 
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looking result was that of Tonry et al. (2003) and Barris et al. (2004) for 
the HZT: the deviation of apparent magnitude from the Iow-CIm, zero-A 
model for several new SNe la at z ~ 1 is roughly the same as that at 
z ~ 0.5, in agreement with expectations based on the results of Riess et 
al. (2001). 

Inspired by his promising results from SN 1997ff, as well as by the 
discovery and study of two distant SNe la (Blakeslee et al. 2003) with 
the HST Advanced Camera for Surveys (ACS), Adam Riess launched an 
extensive HST campaign (the "Hubble Higher- z Supernova Search") to 
closely monitor ~ 6 high-redshift (z ^ 1) SNe la, with the goal of more 
clearly detecting the early epoch of deceleration. A major obstacle was 
the vast amount of HST time required to both discover and obtain follow- 
up photometry of the faint SNe, but this was overcome by arranging to 
"piggyback" on the Great Observatories Origins Deep Survey (GOODS), 
in which ~ 400 HST orbits were being used to obtain deep, multicolor 
ACS images of distant galaxies. The GOODS team obtained their data 
in several distinct epochs with separations of 45 days, and with suitable 
filters, allowing the discovery of high-redshift SNe (e.g., Giavalisco et 
al. 2002; Riess et al. 2004a), generally on the rise. Results from Jha 
(2002) suggested that key light-curve-shape parameters, such as the time 
of maximum and decline rate, could be determined from rest-frame U- 
band light curves, corresponding to observed-frame ACS z-band images, 
allowing for fewer epochs of costly infrared (NICMOS) observations to 
measure the SN luminosity and color. This strategy allowed the requisite 
follow-up images (e.g., Figure 13) and ACS grism spectra to be obtained 
with a relatively modest award of additional HST time. 




Figure 13: SN 2002hp, a high-redshift supernova from the GOODS pro- 
gram using HST. One can see it brightening and subsequently fading 
with time. The assumed host galaxy is at the top of each frame. 

The program was very successful. Riess et al. (2004b; see also Strolger 
et al. 2004) convincingly show that SNe la at z ^ 1 are brighter than 
predicted in simple models invoking systematic effects (Fig. 14). Thus, 
many SN la evolution and dust models are excluded as alternatives to 
acceleration for the observed faintness of z « 0.5 SNe la. The data 
demonstrate that the early universe was indeed apparently decelerating, 
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as expected if A is a relatively recent effect (Fig. 15). The expansion 
of the Universe made a transition between deceleration and acceleration 
(Turner & Riess 2002) at z ~ 0.5, consistent with the new "concordance 
cosmology" of Q,m ~ 0.27, ~ 0.73 from the best CMB measurements 
(Spergel et al. 2003) and (independently) the complete set of available 
SN la data (Riess et al. 2004b). Finally, Riess et al. (2004b) show that 
the value of the dark energy equation-of-state parameter is consistent 
with w = — 1, and that changes in w could not have been very large over 
the past ~ 10 billion years (i.e., dw/dz = w' 0, as expected with A). 
Even better agreement with the A-model expectations was obtained by 
Wang & Tegmark (2004), who analyzed the data in a slightly different 
manner. Thus, although we cannot exclude the possibility that the 
Universe will recollapse in the future, such a "big crunch" (or "gnaB 
giB," which is "Big Bang" backwards) is unlikely to occur in fewer than 
15-30 Gyr (Riess et al. 2004b) if we adopt the linear potential field of 
Kallosh & Linde (2003). 




Figure 14-' SN la residual Hubble diagram comparing cosmological mod- 
els and scenarios for astrophysical dimming (Riess et al. 2004b). Up- 
per panel: SNe la from ground-based discoveries are shown as dia- 
monds, ffiJT-discovered SNe la are shown as filled symbols. Lower panel: 
Weighted averages in fixed redshift bins, which are given only for illustra- 
tive purposes. Data and models are shown relative to an empty universe 
model (O = 0). 
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Figure 15: Kinematic SN la residual Hubble diagram (Riess et al. 2004b). 
Upper panel: Symbols as in Figure 14. Lower panel: Weighted averages, 
as in Fig. 14. Data and kinematic models of the expansion history are 
shown relative to an eternally coasting model, q(z) = 0. Models rep- 
resenting specific kinematic scenarios (e.g., "constant acceleration") are 
illustrated. 



Of course, it is possible to find specific dust or evolution models that 
are not ruled out with our SN la data. For example, a "replenishing 
dust" model represents a constant density of dust that is continually 
replenished at precisely the rate at which it is diluted by the expand- 
ing universe (Goobar, Bergstrom, &; Mortsell 2002; Riess et al. 2004b); 
one could perhaps call this the "steady-state dust model." Since the 
dimming is directly proportional to the distance light traveled and is 
thus mathematically quite similar to the effects of a cosmological con- 
stant, we cannot discriminate this model from the A-dominated model 
in the magnitude-redshift plane. However, the fine-tuning required of 
this dust's opacity, replenishing rate, and velocity ( ^ 1000 km s _1 for 
it to fill space uniformly without adding detectable dispersion) makes it 
unattractive as a simpler alternative to a cosmological constant. More- 
over, the density of intergalactic dust needed to explain the observed 
faintness of SNe la at z ~ 0.5 exceeds the upper limit determined by 
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Paerels et al. (2002) from the absence of a detectable X-ray scattering 
halo around a single distant quasar. 

Another model with this behavior is evolution that is proportional to 
look-back time (Wright 2001). While possible, such dimming behavior, 
especially if in the form of luminosity evolution, seems implausible. We 
may expect evolution (or dust production) to be coupled to the observed 
evolution of stellar populations, galaxies morphologies, sizes, large-scale 
structure, or even chemical enrichment. None of these known varieties of 
evolution are largely completed by z = 0.5 starting from their properties 
at z = 0; quite the contrary, most of them have hardly begun, looking 
back to z = 0.5. As mentioned previously, a strong empirical argument 
against recent luminosity evolution is the independence of SN la distance 
measurements on gross host morphology (Riess et al. 1998b; Sullivan et 
al. 2003). The range of progenitor formation environments spanned by 
SNe la in early-type and late-type hosts greatly exceeds any evolution 
in the mean host properties between z = and z = 0.5. In the end, 
however, the only "proof" against astrophysical contamination of the 
cosmological signal from SNe la is to test the results against those of 
other observations, independent of SNe la. 

4.5 Type la Supernova Rates and Progenitors 

The rate of SNe la as a function of redshift is important to know for 
a number of reasons. For example, it affects the derived chemical evolu- 
tion of galaxies (e.g., Maoz & Gal- Yam 2004; Matteucci & Recchi 2001, 
and references therein), and it can be used to set constraints on the star 
formation history of the Universe (Gal- Yam & Maoz 2004). Moreover, 
the distribution of "delay times" between the formation of the progenitor 
star and the explosion of the supernova can be determined from an accu- 
rate census of the SN la rate at different cosmic times, thereby providing 
clues to the physical nature of the progenitors of SNe la (Madau, Delia 
Valle, & Panagia 1998; Dahlen & Fransson 1999; Yungelson & Livio 
2000). There have been many measurements of the rate of very nearby 
SNe la (Cappellaro, Evans, & Turatto 1999, and references therein), 
and at z « 0.1 several studies have been conducted (Hardin et al. 2000; 
Strolger 2003; Reiss 2000). These form the basis of comparison for the 
measurements at higher redshifts. 

From observations of three SNe la at z ~ 0.4, Pain et al. (1996) 
estimated a rate of SNe la equivalent to 34 yr _1 deg -2 , with a la un- 
certainty of a factor of ~ 2, for objects found in the range of 21.3 < 
R < 22.3 mag. A more recent estimate by Pain et al. (2002a) is 
(1.53 ± 0.3) x 10~ 4 h 3 Mpc -3 yr" 1 at a mean redshift of 0.55 (where 
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h = Ho/(WO km s 1 Mpc 1 ). Cappellaro et al. (1999) report a nearby 
SN la rate of 0.36 ± 0.11 h 2 SNu (where 1 SNu = 10 10 L B& per century), 
or (0.79±0.24) x 10 -4 h 3 Mpc -3 yr" 1 , using a quite uncertain local lumi- 
nosity density of p = 2.2 x 10 8 h L Bq Mpc" 3 . Pain et al. (2002a) claim 
to see a very modest increase in the rate of SNe la with redshift, perhaps 
tracking the star formation rate which Wilson et al. (2002) estimate as 
being proportional to (1 + z) 1 ' 7 in a flat universe with Qm = 0.3. 

Based on a new sample of SNe la observed during the 1999 HZT 
campaign, an independent estimate of the SN la rate at a mean redshift 
of 0.46 was made by Tonry et al. (2003): (1.4 ± 0.5) x 10" 4 h 3 Mpc" 3 
yr" 1 . This is in excellent agreement with the results of Pain et al. 
(2002a), and it is not inconsistent with the local rates from Cappellaro 
et al. (1999), particularly given the uncertainty in the local luminosity 
density. On the other hand, from a different sample of SNe la (see 
below), Dahlen et al. (2004) find a somewhat higher rate of SNe la 
at z ~ 0.4 and suggest that the previously published rates may have 
been underestimated due to systematic effects. Indeed, B. Barris (2004, 
private communication) has found evidence that the true rate of SNe la 
in the 2001 HZT campaign is higher than that reported by Barris et al. 
(2004); similar biases may have affected the 1999 HZT campaign and 
other previous searches. 

Although Tonry et al. (2003) were incomplete in their counts of SNe la 
at z ~ 1, they do not believe that the rate of SNe la closely tracks the 
star formation rate. Application of the star formation rate from Wilson 
et al. (2002) would suggest that the SN la rate at z ~ 1.1 should be 
three times as great as occurs locally and nearly twice as great as the 
rate at z = 0.46. In this case, Tonry et al. (2003) should have discovered 
16 SNe la deg" 2 in their search, yet they only found 4 SNe la deg -2 . 
It is their impression that the constant rate per volume is closer to the 
truth. 

The more recent results of Dahlen et al. (2004) conflict with this 
conclusion. They determine the rate of SNe la as a function of redshift 
by using data from the Hubble Higher- z Supernova Search of Riess et 
al. (2004b), with a redshift range of 0.2 < z < 1.6. The resulting SN la 
rate at z ~ 1 is a factor of 3-5 higher than previous estimates made at 
lower redshifts (z < 0.5), presumably because the star formation rate at 
z > 1 was substantially higher than that in the past few Gyr. Moreover, 
their data suggest that at even higher redshifts, z > 1, the rate of SNe la 
begins to decrease. They find that the delay time (from progenitor star 
formation to SN la explosion) is likely to be substantial, £ 2 Gyr. In 
addition, assuming a Salpeter (1955) initial mass function and a SN la 
progenitor main-sequence mass range of 3-8 M , it appears that 5-7% of 



32 



the white dwarfs formed from these progenitor stars eventually become 
SNe la (but fewer than 1% of all white dwarfs eventually become SNe la). 

Strolger et al. (2004) use the above results, together with additional 
modeling, to set more quantitative constraints on the nature of the pro- 
genitors of SNe la. Specifically, they use a Bayesian maximum likeli- 
hood test to determine the most likely range of delay times that best 
reproduces the observed redshift distribution of SNe la. They find that 
models requiring a large fraction of "prompt" (less than 2 Gyr) SNe la 
poorly reproduce the observed redshift distribution and are rejected at 
> 99% confidence. Thus, SNe la cannot generally be prompt events, 
nor can they be expected to closely follow the star formation rate his- 
tory Instead, Gaussian models best fit the observed data for mean 
delay times in the range of 3-4 Gyr. This may be most consistent 
with single-degenerate systems in which the white dwarf accretes from a 
main-sequence companion or from a somewhat evolved companion (Livio 
2001), although certain types of double-degenerate models are not yet 
eliminated. Tests conducted by Gal- Yam & Maoz (2004) also conclude 
that the characteristic delay times of SNe la should be large (>l-2 Gyr) 
for similar assumed models of the star formation rate history, but the 
results are not as definitive as those of Strolger et al. (2004) because 
they are based on SN la rates derived from more limited SCP data (Pain 
et al. 2002a). 

4.6 Measuring the Dark Energy 
Equation-of-State Parameter 

Every energy component in the Universe can be parameterized by 
the way its density varies as the Universe expands (scale factor a), with 
p oc a^ l+w \ and w is the component's equation-of-state parameter, 
w = P/(pc 2 ), where P is the pressure exerted by the component. So for 
matter, w = 0, while an energy component that does not vary with scale 
factor has w = — 1, as in the cosmological constant A. Quintessence 
models have w ^ — 1, and generally dw/dz ^ 0). Some really strange 
energies may have w < — 1: their density increases with time (Carroll, 
Hoffman, &; Trodden 2003), leading to a "Big Rip" in which progressively 
smaller bound systems get torn apart! [Riess et al. (2004b) and Caldwell 
et al. (2003) estimate that such a fate will not occur sooner than ~ 20 
Gyr from now, if ever.] Clearly, a good estimate of w becomes the key 
to differentiating between models. 

The CMB observations imply that the geometry of the universe is close 
to flat, so the energy density of the dark component is simply related to 
the matter density by £l x = 1 — Qm ■ This allows the luminosity distance 
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as a function of redshift to be written as 
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c(i + z) r z [i + n g ((i + z) 
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showing that the dark energy density and equation of state directly 
influence the apparent brightness of standard candles. As demonstrated 
graphically in Figure 16, SNe la observed over a wide range of redshifts 
can constrain the dark energy parameters to a cosmologically interesting 
accuracy. 

But there are two major problems with using SNe la to measure w. 
First, systematic uncertainties in SN la peak luminosity limit how well 
Dl(z) can be measured. While statistical uncertainty can be arbitrarily 
reduced by finding thousands of SNe la, intrinsic SN properties such as 
evolution and progenitor metallicity, and observational limits like pho- 
tometric calibrations and K-corrections, create a systematic floor that 
cannot be decreased by sheer force of numbers. We expect that system- 
atics can be controlled to at best 3%, with considerable effort. 

Second, SNe at z > 1.0 are very hard to discover and study from the 
ground. As discussed above, both the HZT and the SCP have found 
a few SNe la at z > 1.0, but the numbers and quality of these light 
curves are insufficient for a w measurement. Large numbers of SNe la 
at z > 1.0 are best left to a wide-field optical/infrared imager in space, 
such as the proposed Supernova/ Acceleration Probe (SNAP; Nugent et 
al. 2001) satellite. 

Fortunately, an interesting measurement of w can be made at present. 
The current values of Qm from many methods (most recently WMAP: 
0.27; Spergel et al. 2003) make an excellent substitute for those ex- 
pensive SNe at z > 1.0. Figure 16 shows that a SN la sample with a 
maximum redshift of z = 0.8, combined with the current 10% error on 
£Im, will do as well as a SN la sample at much higher redshifts. Within 
a few years, the Sloan Digital Sky Survey and WMAP will solidify the 
estimate of VLm and sharpen w further. 

Both the SCP and the HZT are involved in multi-year programs to 
discover and monitor hundreds of SNe la for the purpose of measuring 
w. For example, the HZT's project, ESSENCE (Equation of State: Su- 
pErNovae trace Cosmic Expansion), is designed to discover 200 SNe la 
evenly distributed in the 0.2 < z < 0.8 range (Smith et al. 2002; Gar- 
navich et al. 2002; http://www.ctio.noao.edu/wproject). The CTIO 
4-m telescope and mosaic camera are being used to find and follow the 
SNe by imaging on every other dark night for several consecutive months 
of the year. Keck and other large telescopes are being used to get the 
SN spectra and redshifts. Project ESSENCE will eventually provide an 
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estimate of w to an accuracy of ~10% (Figure 17). Even larger num- 
bers of high-redshift SNe la are being found during the ongoing CFHT 
Legacy Survey (http://www.cfht.hawaii.edu/Science/CFHLS; Pain et 
al. 2002b), providing an independent sample with which to measure 
the value of w. Within the next few years, Pan-STARRS (Kaiser et al. 
2002; http://poi.ifa.hawaii.edu) should discover and follow thousands of 
SNe la. 




[)„=1.0-0„ 

Figure 16 (left): Constraints on Q x and w from SN data sets collected 
at z = 0.2 (solid lines), z = 0.7 (dashed lines), and z = 1.6 (dash-dot 
lines). The shaded area indicates how an independent estimate of Q m 
with a 10% error can help constrain w. 

Figure 17 (right): Expected constraints on w with the desired final 
ESSENCE data set of 200 SNe la, 30 of which (in the redshift range 
0.6 < z < 0.8) are to be observed with HST. The thin lines are for SNe 
alone while the thick lines assume an uncertainty in Ojvf of 7%. The 
final ESSENCE data will constrain the value of w to ~10%. 

Farther in the future, the plethora of SNe la to be found and studied 
by the proposed SNAP satellite (Nugent et al. 2001), the NASA/DOE 
Joint Dark Energy Mission (JDEM), the Large-area Synoptic Survey 
Telescope (the "Dark Matter Telescope"; Tyson & Angel 2001), and 
similar large-scale projects could reveal whether the value of w depends 
on redshift, and hence should give additional constraints on the nature of 
the dark energy. High-redshift surveys of galaxies such as DEEP2 (Davis 
et al. 2001), as well as space-based missions to map the CMB (Planck), 
should provide additional evidence for (or against) A. Observational 
cosmology promises to remain exciting for quite some time! 
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